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p ■ ABSTRACT 

I \ HD 133880 is a rapidly rotating Bp star (y sin i ~ 103 km s""'^) and is host to one of the 

O . strongest magnetic fields of any Ap/Bp star. A member of the Upper Centaurus Lupus 

association, it is a star with a well-determined age of 16 Myr. Twelve new spectra, 
^ \ four of which are polarimetric, obtained from the FEROS, ESPaDOnS, and HARPS 

instruments, provide sufficient material from which to re-evaluate the magnetic field 
and obtain a first approximation to the atmospheric abundance distributions of He, 
O, Mg, Si, Ti, Cr, Fe, Ni, Pr, and Nd. 
^ ' An abundance analysis was carried out using ZEEMAN, a programme which syn- 

C/^ \ thesises spectral line profiles for stars with permeating magnetic fields. The magnetic 

field structure was characterised by a colinear multipole expansion from the observed 
variations of the longitudinal and surface fields with rotational phase. Both magnetic 
hemispheres are clearly visible during the stellar rotation, and thus a three-ring abun- 
CO . dance distribution model encompassing both magnetic poles and magnetic equator 

^^ ' with equal spans in co-latitude was adopted. 

Using the new magnetic field measurements and optical photometry together with 
previously published data, we refine the period of HD 133880 to P = 0.877476 ± 
0.000009 days. Our simple axisymmetric magnetic field model is based on a pre- 
dominantly quadrupolar component that roughly describes the field variations. Using 
^ . spectrum synthesis, we derived mean abundances for O, Mg, Si, Ti, Cr, Fe, and Pr. 

j^ \ All elements, except Mg, are overabundant compared to the Sun. Mg appears to be 

" " approximately uniform over the stellar surface while all other elements are more abun- 

dant in the negative magnetic hemisphere than in the positive magnetic hemisphere. 
In contrast to most Ap/Bp stars which show an underabundance in O, in HD 133880 
this clement is clearly overabundant compared to the solar abundance ratio. 

In studying the Ha and Paschen lines in the optical spectra we could not unam- 
biguously detect information about the magnetosphere of HD 133880. However, radio 
emission data at both 3 and 6 cm suggests that the magnetospheric plasma is held 
in rigid rotation with the star by the magnetic field and further supported against 
collapse by the rapid rotation. Subtle differences in the shapes of the optically thick 
radio light curves at 3 and 6 cm suggest that the large-scale magnetic field is not fully 
axisymmetric at large distances from the star. 
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1 INTRODUCTION 

HD 133880 (=HR 5624) is a rapidly rotating late B-type 
chemically peculiar (Bp) star with the Si A4200 pecuUar- 
ity. It exhibits an unusual magnetic field. The mean line- 
of-sight magnetic field {B^) is very s trong an d observed to 
vary firom about -4 to +2 kG (Lands treetlll990l ). Unlike most 
chemically peculiar stars, HD 133880 has a field that is pre- 
dominantly quadrupolar as opposed to dipolar. HD 133880 
is a photometric variable with variations on the order of 
0.15 mag in the t/-band, which m ay be the result of the large 
magnetic field (| Waelkend 1 1985^ . The magnetic field serves 
to create a "patchy" distribution of elements on the surface 
of the star. This non-uniform abundance distribution, and 
its associated line blocking and backwarming, may explain 
the photometric variations (Rvabchikova 1991 1 . Specifically, 
several elements are distributed non-uniformly over the stel- 
lar surface in a non-axisymmetric pattern about the rota- 
tion axis. As the star rotates, the observed magnetic field 
strength (Bz) varies and over/underabundances of elements 
are detected. The most striking anomalous abundances in 
Bp stars are often found for Cr and rare-earths which can 
be as much as 10^ and 10"^"^ ove rabundant compared to the 
Sun, respectively (jRvabchikoval ll991V The resultant vari- 
ability is explained by the oblique rotator model: the rota- 
tion and magnetic field axes are at angles i and j3 to the 
line -of-sight and rotat ion axis, respectively. 

iLandstreetl ([1990|) used longitudinal magnetic field mea- 
surements obtained fro m H/3 in c onjunction with photo- 
metric data reported by IWaelkeng ([1985) to deduce a pe- 
riod of P = 0.877485 ± 0.00002 days for HD 133880. 
iLim. Drake fc Linskvl (|1996l ) did a preliminary analysis of 
the rotational modulation of 3 cm and 6 cm radio emission 
from HD 133880. The radio emission variation had both 
broad and narrow peaks that correlated with the maximum 
extrema of the dipolar, Bd, and quadru polar, Bg, contribu - 
tions to the magnetic field reported bv ILandstreetl l|l99(J l. 
resp ectively. 

iLandstreet et all (|2007l ') confirmed HD 133880 as a 
member of the Upper Cen Lup association and thus it is 
a star with a well known age of logt = 7.20 ± 0.10 (yr). 
HD 133880 is a very young star, havin g completed only 
about 5% o f its main sequence lifetime (JLandstreet et al.l 
I2OO7I . 120081 ). Their preliminary analysis of many physi- 
cal characteristics concluded that T^s ~ 12000 ± 500 K, 
log L/Lq = 2.10 ± 0.1, and M/Mq = 3.20 ± 0.15. (Our 
analysis (see below) has shown that the Tee is likely closer 
to 13000 K). Table 1 summarises the physical properties of 
HD 133880. 

There are few hot chemically peculiar stars with broad 
spectral features for which abundance analyses have been 
carried out. Because of the very complex, strong and un- 
usual magnetic field, HD 133880 is a worthwhile candidate to 
study in detail. It is also interesting from the fact that it is a 
star with a well known age. HD 133880 is one in a large sam- 
ple of clust er Ap/Bp stars first co r npiled by iBamulo et al.l 
l|2006t) and ILandstreet et al.l (|2007l . I2OO8I ') to study the time 



Hawaii Telescope (CFHT) which is operated by the National Re- 
search Council of Canada, the Institut National des Science de 
I'Univers of the Centre National de la Recherche Scientifique of 
France, and the University of Hawaii. 



Table 1. Summary of stellar, wind, magnetic and magnetospheric 
properties of HD 133880. 



Spectral type 


B8IVp Si A4200 


SIMBAD & BSC 


Teff (K) 


13000 ± 600 


This paper 


log g (cgs) 


4.34 ± 0.16 


This paper 


R. (Ro) 


2.01 ± 0.32 


This paper 


vsini (kms^-"^) 


103 ± 10 


This paper 


P(d) 


0.877476 ± 0.000009 


This paper 


log(L*/L0) 


2.02 ± 0.10 


This paper 


M* (Mq) 


3.20 ± 0.15 


Landstreet et al. (2007) 


logt (Myr) 


7.20 ± 0.10 


Landstreet et al. (2007) 


Bd(G) 


-9600 ± 1000 


This paper 


Sq(G) 


-23000 ± 1000 


This paper 


Boct (G) 


1900 ± 1000 


This paper 


i{°) 


55 ± 10 


This paper 


/3(°) 


78 ± 10 


This paper 


r]. 


^10^ 


This paper 


W 


0.3 


This paper 


'''spin 


11 Myr 


This paper 



evolution of magnetic fields in Ap stars. The current study 
is part of an effort to understand empirically how chemical 
abundances evolve with time and s tellar age in the ma gnetic 
chemically peculiar stars. Recentlv lBailev et al.l (|201ll l char- 
acterised one star in this sample, the Ap star HD 318107 
(=NGC 6405 77). HD 133880 is a second star of great inter- 
est from that sample. 

This paper will discuss the efforts to characterise the 
magnetic field and abundance distributions of several ele- 
ments for HD 133880. The following section discusses the 
polarimetric and spectroscopic observations used as well as 
the deduced longitudinal fields; Sect. 3 discusses improve- 
ment of the rotational period ; Sect. 4 describes the determi- 
nation of physical parameters; Sect. 5 addresses the observed 
line profile variations; Sect. 6 & 7 discuss the magnetic field 
and abundance models; Sect. 8 reports the abundances ob- 
tained; Sect. 9 discusses the Ha and radio magnetosphere; 
and Sect. 10 summarises the presented work. 



2 OBSERVATIONS 

2.1 Polarised spectra 

High-resolution spectropolarimetric (Stokes I and V) ob- 
servations of HD 133880 were collected with ESPaDOnS at 
the 3.6m Canada-France-Hawaii Telescope and HARPSpol 
on the ESQ 3.6m telescope. ESPaDOnS and HARPSpol 
are both high resolution fibre-fed spectropolarimeters. Three 
ESPaDOnS observations were obtained on August 3, 2010, 
and on July 14 and 15, 2011 (the latter in the context of the 
Magnetism in Massive Stars Large Programme (MiMeS)), 
and one HARPSpol spectrum was obtained on May 3, 2010. 
Each ESPaDOnS spectropolarimetric sequence consisted of 
four individual subexposures taken in different retarder con- 
figurations, whereas the HARPSpol spectrum results from 
the combination of four independent spectra, each comprised 
of four subexposures, obtained over the course of about 9 
minutes. The ESPaDOnS measurements cover a spectral 
range from 370 nm to 1050 nm at a spectral resolution of 
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Table 3. Log of unpolarised spectra. Listed are the instrument 
used, the Hehocentric Julian date, the peak SNR per 1.8 kms~ 
velocity bin, the exposure time and the phase computed using the 
ephemeris of Eq. (2). 



Instrument 


HJD 


SNR 


Ccxp (s) 


Phase 


ESPaDOnS 


2455412.727 


180 


386 


0.774 


ESPaDOnS 


2455414.728 


180 


394 


0.054 


FEROS 


2454867.846 


130 


143 


0.809 


FEROS 


2454868.818 


130 


103 


0.917 


FEROS 


2455031.608 


250 


122 


0.438 


FEROS 


2455599.820 


240 


141 


0.990 


FEROS 


2455600.828 


240 


151 


0.139 


FEROS 


2455601.868 


300 


121 


0.325 



R~65 000, while the HARPSpol observation covers 370 nm 
to 690 nm at a spectral resolution of R~100 000. 

From each set of four subexposures we derived Stokes / 
and Stokes V spectra following the double-ratio procedure 
described bv lDonati et al.l (jlOOTj), ensuring in particular that 
all spurious signatures are removed at first order. Null po- 
larization spectra (labeled A'^) were calculated by combining 
the four subexposures in such a way that polarisation can- 
cels out, allowing us to verify that no spurious signals are 
present in the data (see Donati et al. 1997 for more details on 
the definition of N). All ESPaDOnS frames were processed 
using the CFHT's Upena pipeline, feeding the auto mated 
reduction package Libre ESpRIT (|Donati et aLlliggTl ). The 
HARPSpol spectrum was r educed using the REDUCE code 
l|Piskunov fc Valentil 12003 ) . The details of this reduction 
can be found in recen t paper s by iMakaganiuk et al.l ()2011h 
and iKochukhov et al.l (|2011bl ). The peak signal-to-noise ra- 
tios (SNRs) per 1.8 kms~^ velocity bin in the reduced ES- 
PaDOnS and HARPSpol spectra range from 449-697. 

The log of spectropolarimetric observations is presented 
in Table m 



2.2 Unpolarised spectra 

High resolution spectroscopic observations (Stokes I) were 
obtained from ESPaDOnS at CFHT and ESO's FEROS in- 
strument at the 2.2m telescope located at La Silla. FEROS is 
an Echelle spectrograph that has a wavelength coverage from 
350 nm to 920 nm with a spectral resolution of R~48000. 
FEROS is fed by two fibres (object and sky fibres), which 
in principle can be used to subtract the contribution of the 
sky from t he object spec trum. A more detailed description 
is given bv lKaufer et all (|l999l ). The two ESPaDOnS obser- 
vations were taken on August 4 and 6, 2010, and the six 
FEROS spectra were collected on February 5, 6 and July 
19, 2009 as well as February 7, 8 and 9, 2011. 

The ESPaDOnS data were reduced using Libre-Esprit 
in a manner similar to that described above. The FEROS 
data were reduced using the data-reductio n software (DRS) 
implemented under MIDAS (JKaufer et al.l 1999). The log of 
unpolarised spectra is presented in Table [3| 



2.3 Longitudinal Magnetic Field Measurements 

The circular polarisation observations were analyzed using 
the multiline analysis technique Least Squares Deconvolu- 
tion (LSD; Donati et al. 1997). LSD combines the informa- 
tion from essentially all metallic and He lines in the spec- 
trum by means of the assumption that the spectrum can 
be reproduced by the convolution of a single "mean" line 
profile (the LSD profile) with an underlying spectrum of 
unbroadened atomic lines of specified line depth, Lande fac- 
tor, and wavelength (the line mask, as described by Wade 
et al., 2000). This process allows the computation of single 
averaged Stokes I and V profiles with much higher signal- 
to-noise ratios than those of the individual lines, dramati- 
cally improving the detectability of Zeeman signatures due 
to stellar magnetic fields. 

The atomic data were take n from the Vienna 
Atomic Line Database (VALD) llPiskun ov ct al.' Il995l : 
iRvabchikova et al.l Il997l : iKupka et all |l999, 2000). The 
VALD line list used in this analysis is characterised by the 
stellar effective temperature T^s and surface gravity log g re- 
ported in Table 1, abundances as derived from the detailed 
abundance analysis (see Sect. 8), and an unbroadened line 
depth threshold equal to 1% of the continuum intensity /c, 
yielding approximately 10000 lines in the ESPaDOnS spec- 
tral window and 8000 lines in the HARPSpol window. The 
line depths were then interactively adjusted so that the line 
depths predicted by the LSD convolution model matched 
the observed depth of the spectral absorption lines as well 
as possible, thus improving the fit between the spectrum and 
the LSD model. Because Stokes V profiles are detected in 
many individual spectral lines of HD 133880, we verified (by 
examining simultaneously the Stokes V and / profiles) that 
this adjustment procedure also resulted in a better agree- 
ment between the observed and computed Stokes V spec- 
trum. 

The mean longitudinal magnetic field was measured by 
computing the first-order moment of the Stokes V LSD pro- 
file within the line according to the expression: 



(B^) = -2.14 X 10^ 



^vV{v)dv 



>^gcscj [Ic - I{v)]dv' 



(1) 



where (Bz) is in G, goff is the effective Lande factor, A is 
the mean wavelength of the LSD line in A, and v is the 
velocity within the profile measured relative to the centre 
of gravity (Mathys et al., 1989; Donati et al. 1997; Wade 
et al., 2000). The uncertainties associated with {B^} were 
determined by propagating the formal (photon statistical) 
uncertainties of each pixel through Eq. (1). The integration 
ranges employed in the evaluation of Eq. (1) associated with 
each Stokes profile were selected individually through visual 
inspection so as to include the entire span of the (highly 
variable) Stokes / profile. 



2.4 ATCA radio continuum measurements 

HD 133880 was observed with the Australia Telescope Com- 
pact Array (ATCA) at 6 cm and 3.5 cm wavelength simul- 
taneously on February 12, 14 and 16, 1995. On each day, the 
observation spanned ~10 hours, during which time we ob- 
served HD 133880 continuously apart from short (~3 mins) 
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Table 2. Log of ESPaDOnS and HARPSpol observations of HD 133880. Listed are the instrument used, the Heliocentric Julian date of the 
midpoint of the observation, total exposure time, the peak signal-to-noise ratio per 1.8 kms~ velocity bin, the phase of the observation 
(according to Eq. 2), the evaluation of the detection level of a Stokes V Zeeman signature (DD=definite detection, MD=marginal 
detection, ND=no detection), and the derived longitudinal field and longitudinal field detection significance z from both V and N. In 
no case is any marginal or definite detection obtained in the N profiles. 















V 




N 




Instrument 


HJD 


Ccxp 


SNR 
pix~^ 


Phase 


Detect 


Bi it ag 
(G) 


z 


Bi ± (tb 
(G) 


z 


ESPaDOnS 


2455411.732 


720 


629 


0.640 


DD 


2029 ± 38 


53.4 


-51 ±31 


1.6 


ESPaDOnS 


2455756.795 


720 


697 


0.878 


DD 


-2103 ± 38 


55.3 


43 ±23 


1.9 


ESPaDOnS 


2455757.877 


720 


573 


0.111 


DD 


-3671 ± 43 


85.4 


-23 ± 26 


0.88 


HARPSpol 


2455319.726 


520 


449 


0.787 


DD 


650 ± 62 


10.5 


-15 ±58 


0.26 



scans of a secondary calibrator at regular (^^20 mins) in- 
tervals. In this way, we were able to attain full coverage 
in rotation phase for the star, with repeated coverage (on 
two separate days) for just over half a rotation phase to 
check where any observed variations were repeatable (as in- 
deed t hey were, as shown in Fig. 2 of iLim. Drake fc Linskvl 
l| 19961 )). The primary (absolute flux density) calibrator used 
was PKS1934-638, and the secondary (complex gain) cali- 
brator PKS1458-391. We reduced the data in the standard 
fashion using MIRIAD, being careful to omit data with poor 
phase coherence (typically at the beginning and end of each 
day when the source was low in the sky). 



3 ROTATION PERIOD 



A per iod of 0.87746 ± 0.00001 d was reported bv lWaelken^ 
l|l985h based on Geneva photometric observations. Land- 
street (1990) recommended an adjustment of the period to 
0.877485 ± 0.000020 d to b ring the ma gnetic measurements 
reported by Borr a fc Landstreetl (|l97a) into agreement with 
those that he had acquired in 1987 and 1988. 

When we phase the accumulated magnetic mea- 
surements of HD 133880 _(i.e. the three measurements 



of Borra fc Landstroetj (|l975l ). the 12 measurements of 



iLandstreet (1990 ). and the 4 new measurements presented 
here) according to Landstreet's (1990) ephemeris, we ob- 
serve a systematic offset of ~ 0.1 cycles between the new 
measurements and those published previously. Such an off- 
set could be interpreted as the consequence of a small error 
in the assumed period. However, given the presence of im- 
portant abundance non-uniformities on the surface of HD 
133880, it is likely that the shape, and possibly the ampli- 
tude of the longitudinal field variation derived from metallic 
lines (i.e. our new LSD measurements) and those derived 
from hydrogen lines (i.e. the previously published measure- 
ments) may differ. Given the sparse phase coverage of our 
new measurements, it is therefore ambiguous whether this 
apparent offset is truly a consequence of a period error. 

In an attempt to resolve this uncertainty, we exam- 
ined the photometric measurements obtained by the Hip- 
parcos mission (ESA, 1997). One hundred iifp-band mea- 
surements of HD 133880 are reported in the 1997 edition 



of the cataloguelj The periodogram of those measurements 
in the range 0.5-1.5 days shows a clear, unique peak at 
0.877479±0.000030 d (The uncertainty on the period, which 
corresponds to la confidence, is computed assuming Gaus- 
sian statistics by calculating the variation of the reduced y^ 
of a sinusoidal fit to the phased measurem ents in the vicin- 
ity of the best-fit period (e.g. fress, Flann erv fc Teukolskvl 
1 19861 )). A similar analysis applied to Waelkens' (1985) 37 V- 
band measurements yields 0.877459J;o:oooo6 d. Both of these 
periods are consistent within the (Icr) error bars. To improve 
the precision of the period, we combined the two datas ets. 
First, we used the transformation reported by iHarmaned 
1 19981 ) to convert the Hp measurements to Johnson V-band. 
We then performed the same period analysis on the com- 
bined data set (spanning nearly 15 years), obtaining a period 
of 0.8774731 ± 0.0000060 d. 

Therefore the photometric data are indicati ve of a pe- 
riod o f 0.87 74731 d, intermed iate between those of lWaelkenj 
1 19851 ) and ILandstreet l|l990l ). 

We note that if we phase the magnetic data with 
the above photometric period, a slightly longer period of 
0.877476 d is necessary to bring the new magnetic measure- 
ments into better agreement with the older measurements. 
Finally, this period also removes the p hase off set between the 
observed radio fiux variation noted bv lLim. D rake fc Linskvl 
(| 19961 ) and the magnetic extrema. 

Based on these results, we adopt the following 
ephemeris for HD 133880: 



jD/fQ" 



2445472.000(10) +0.877476(9) • E, 



(2) 



where we have adopted a zero point corresponding to min- 
imum photometric brightness (which occurs simultaneously 
with the minimum of the longitudinal magnetic field). We 
determined this (3(t) uncertainty by examining the x^ of the 
fit to the photometric data about this best-fit value. 



vizier . u-strasbg . fr/viz-bin/VizieR-S?HIPy.2074066) 



HD 133880 5 




Figure 1. Phased V-band llWaelkenj lll985l) : filled squares) and 
i?p-band (ESA 1997; open diamonds) photometry. The Hipparcos 
photometry has been scaled to the Johnson V system using the 
transformation of lHarmaned 1119981 ). 



4 DETERMINATION OF PHYSICAL 
PARAMETERS 

4.1 Effective Temperature 

The enhanced metal abundances and the effects of 
the magnetic field have to be taken into account 
when determining the effecti v e temperatures of Ap stars 
(jStepien fc DominiczakI 1 19891 : iHauck fc Kunzlil 1 19961 ). As 
a result, effective temperatures determined using Geneva 
and/or Stromgren uvbyP photometry using normal star cal- 
ibrations must be corrected to t he Ap stars temperature 
scale. However, iHauck fc NorthI l| 19931 ) suggest that stars 
that are either He-weak or He-strong do not require that 
their temperatures be corrected. 

HD 133880 has available both Geneva and Stromgren 
uvbyP photometry for the determination of effective 
temperature. We have used both in this analysis. 
For the Geneva photometry we have used th e FOR- 
TRAN programme desc ribed bv iKunzli et all ([l993), 
and discussed by iLandstreet et al.l ( 20071 ). to obtain 
an estimate of the effective temperature. For the 
Stromgren uvbyP photometry we used the FORTRAN code 
"UVBYBETANEW-a p" which is a modified version of 
"UV B YBETANEW" (iNapiwotzki. Schoenberner fc Wenskj 



1 19931 : iMoon fc Dworetskvl Il985h that corrects the effec- 
tive temperature to the Ap temperature scale. However, 
in this analysis no correction was applied to the temper- 
ature because HD 13 3880 is classified as a He-weak star 
IjHauck fc Northlll993l ). which is confirmed in this paper (see 
Sect. 8). 

We found Tctt = 13300 K and 12700 K using uvbyP 
and Geneva photometry respectively. The scatter of these 
measurements about their mean is of the order of 300 K. 
Given the intrinsic uncertainties in determining the effec- 
tive temperature for a magnetic Ap star in addition to the 
global uncertainties in the applied methods we optimisti- 
cally adopt an uncertainty of approximately ±600 K. Both 



the Geneva and uvbyP photometries are rotationally av- 
eraged values and their variations are included within the 
quoted uncertainties. We use the average of the computed 
temperatures, adopting an effective temperature of Tcfi = 
13000 ± 600 K for HD 133880. If the Ap correction to the 
temperature is applied, a value of T^h = 12440 K is found, 
whic h is within the the uncertainty of our adopted value. 

INetopil et al.l l|2008h performed an analysis on the tem- 
perature calibration of chemically peculiar stars. Included 
in their analysis was HD 133880. Using photometry, their 
calibrated temperature corrections suggest that the Tctt is 
11930 ± 210 K. This agrees, within uncertainties, to our 
derived Tcff if we corrected our derived temperature to the 
Ap temperature scale . In our analysis , we fa vour the higher 
Tcft based on work bv lHauck fc NorthI ll 19931) and the b etter 
quality fits to H/? using the model fits of lKuruca (|l979l ) with 
the higher of these two temperatures. We do, however, com- 
pare the de rived abundances of all elements using the Tofi 
reported bv lNetopil et al.l (|2008l ) of 12000 K to our adopted 
value of 13000 K (see Sect. 8). 



4.2 Luminosity 

In addition to Tctt and the V magnitude of HD 133880 
(see Table [!}, we also require a bolome tric correction in 
order to calculate the stellar luminosity. ILandstreet et al.l 
((20071) derived a set of bolometric corrections as a function 
of Teff that can be applied to Ap stars. Ap p lying this cor- 
rection (see Equation 1 in ILandstreet et al.l (|2007|)) results 
in log L/Lq = 2.02 ± 0.10. Here, the uncertainty was esti- 
mated in the same ma nner as described in Sect. 1 and by 
ILandstreet et al.l (|2007h . Note that we used E(B - V) = 0.00 
to compute the luminosity, which is confirmed by looking at 
the colour-colour diagram of (U - B) versus (B - V) 



4.3 Other Parameters 

Based on values of Toff and luminosity determined in 
Sect. 4.1 and 4. 2, the stellar radius is found to be i? = 
2.01 ± 0.32 _Rq. ILandstreet et all (|2007l ) report a mass for 
HD 133880 of M = 3.20 ± 0.15 Mq (consistent with the 
values of T^s and luminosity in this paper), which can be 
used with the derived radius to estimate the surface gravity, 
log g. We find log g — 4.34 ± 0.16, consistent with young 
age. 

From our adopted period of P = 0.877476 ± 
0.000009 days (Sect. 3) and «sinj = 103 ± 10 kms"^ (see 
below) we can estimate the inclination angle of the rotation 
axis to the line of sight, i, using the equation 



{v sin i)P 



(3) 



50.QR ' 

where vsini is in kms~^ and R in solar radii, to obtain a 
value of i = 63° ± 18° . All the physical parameters derived 
are summar ised in Table [T] and are in agreement with those 
reported bv ILandstreet et al.l (|2007l ). 



5 LINE PROFILE VARIATIONS 

To characterise the line profile variability of HD 133880 we 
measured the equivalent width variations for a number of 
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Figure 2. Phased equivalent widths measurements for the in- 
dicated spectral lines measured from the ESPaDOnS (black 
squares), FEROS (red triangles) and HARPSpol (blue circles) 
datasets. 



spectral lines. The spectral lines were re-normalised to the 
surrounding continuum regions and the equivalent widths 
were computed by numerically integrating over the line pro- 
file. A single uncertainty value was estimated for each pixel 
from the RMS scatter in the continuum regions surrounding 
the line profile and the resulting uncertainties in the equiv- 
alent width measurements were found by adding the pixel 
uncertainties in quadrature. 

In Fig. [2] we show equivalent width measurements 
phased to the ephemeris given in Eq. ([2]) for select lines of 
titanium, iron, chromium and silicon. The equivalent width 
curves are all similar and show strong, approximately si- 
nusoidal variations consistent with the rotation period. For 
each element, the equivalent width reaches a minimum at 
about phase 0.5. The largest variation is seen in the Cr ll 
4588 A line with the equivalent width measurements al- 
most doubling in value. Similar equivalent width changes 
are found in the Till 4533 A and Fell 4583 A lines, but a 
much smaller amplitude of the variation is measured in the 
Sill 5056 A doublet. Clear variations are also seen in other 
spectral lines of these elements. 

The equivalent width variations for each of the elements 
plotted in Fig. [2] appear to be the result of absorption fea- 
tures travelling through the line profile. The net result is a 
variation in the total depth, as well as shape, of the line, 
as demonstrated by the profiles shown in Fig. [3] The Si 
line plotted in Fig. [3] has the simplest variations, with what 
appears to be a single feature travelling from negative to 
positive velocities. The variations in the the other lines are 
phased similarly, but the line profiles are generally more 
complex. 



6 MAGNETIC FIELD MODEL 



Before any meaningful abundance analysis can be done, 
we need to establish an appropriate magnetic field model. 
iLandstreeJ ((1990) derived a best fit magnetic field model 
based on longitudinal field measurements of H/3 for 
HD 133880. He adopted a magnetic field geometry with the 
inclination angle of the line of sight to the rotation axis, i, 
equal to the angle of the magnetic field axis to the rota- 





4532 4533 4534 4535 4536 4537 
wavelength (angstroms) 



CrII 4588/4592 



4582 4583 4584 4585 4586 
wavelength (angstroms) 




4588 4590 4592 4594 

wavelength (angstroms) 



5053 5054 5055 5056 5057 5058 5059 
wavelength (angstroms) 



Figure 3. Observed metallic line profiles for the indicated lines. 
Shown are the continuum-normalised spectra, displayed in such 
a way that the continuum of each spectrum is plotted at a po- 
sition on the vertical axis that corresponds to the phase of the 
observation. 



tion axis, /3 with a value of 90° {i = P — 90°). The dipo- 
lar and quadrupolar components were Bd — —8125 G, and 
Bq = —10900 G respectively. (Note that we have corrected 
negative sig n errors in both co mponents for the derived field 
reported bv lLandstreetl (|l990l )). From this model, it was es- 
tablished that the predominant magnetic field distribution 
over the stellar surface is quadrupolar with (Bz) variations 
that are clearly not sinusoidal. Further, since a change in the 
sign of the longitudinal field is observed, i+P must be greater 
than 90°. However, Landstreet's adopted geometry is barely 
consistent with our calculated inclina tion of i = 63° ± 18° 
(see Sect. 4). Also, iLandstreed (|l990l ) had no information 
about the variations of the surface magnetic field modu- 
lus, (B), making it difficult to constrain the field further. 
iLandstreeti l|l99d ) pointed out that there is not a unique 
model for the observed magnetic curve of HD 133880 and 
suggested an alternate geometry of i = 80° and P = 55° that 
produces a similar curve (note that the values of i and /3 can 
be interchanged without affecting the magnetic curve) . This 
second geometry is therefore in much better agreement with 
our derived value of i. Unfortunately, iLandstrect (1990) did 
not provide the values of the dipole and quadrupole compo- 
nents for this alternate geometry. We therefore search for a 
new magnetic field model that agrees well with our derived 
value for i of 63° ±18°. 
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Table 4. Lines used to measure {B). The atomic data were gath- 
ered from the VALD database. 





' 1 ' 


1 ' 1 


' 1 ' 


1 1 


2000 


- 


^^T~- 




- 


1000 


- 


/ 




- 





- 


/ 




- 


fl -1000 

■9 


- 


/ 




- 


-2000 




/ 


i 


y 


-3000 


- \ 1 






V J 


-4000 


1 


1 , 1 


1 


1 



0.25 0.5 0.75 1 

Phase 

Figure 4. The top panel depicts the surface magnetic field modu- 
lus variations as measured from Fe II (red upward facing triangles) 
and Cr II (blue downward facing triangles) with the solid (black) 
curve representing the adopted best-fit model. The bottom panel 
shows the (B^) field variations observed for HD 133880 for our 
adop ted magnet i c field model. The red circles are measurements 
from lLandstreeti l|l990f ) using H/3, the blue triangles are LSD mea- 
surements from ESPaDOnS, the purple diamond HARPS and the 
solid (black) curve the adopted magnetic field model. The mini- 
mum of the (Bz) and (B) curves occur at phases 0.037 and 0.537 
respectively. 



6.1 Surface magnetic field modulus measurements 

To further constrain the magnetic field, information about 
the surface magnetic field modulus (B) is required. To obtain 
this, we cannot rely on Zeeman splitting because rotational 
broadening dominates the spectrum, smearing out any ob- 
servable Zeeman components. Although dominated by rota- 
tion, the magnetic field contributes a m easurable fracti on to 
the width of a (Stokes /) spectral line. IPrestonI l|l971f ) out- 
lines a method to obtain mean surface field measurements 
from rotationally broadened spectral lines. By comparing 
the widths of two spectral lines of a given element that 
have a large and small mean Lande factor (ie. lines that 
are strongly effected by the local magnetic field to ones that 
are not) we can obtain a mean i ngful measurement of the 
magnetic field modulus. IPrestonI (|l97ll ) defines a parameter 
K: 

1/2 

K = ( ,,. „:" :,". ., I , (4) 



(X^Z^)^ - (A4^2), 

where wl and ws are the averages of the measured line 
widths, in cm, of lines with a large and small Lande factor 
respectively, A is the wavelength in cm and z the mean Lande 



Element 


Small 

A (A) 


z 
z 


Large 
A (A) 


^ 


Fe II 
Cm 


4508.288 
4284.188 


0.50 
0.52 


4520.224 
4261.913 


1.34 
1.080 



Table 5. (B) measurements from lines of Cr II and Fe II for 
spectra from ESPaDOnS and HARPSpol. Phases are computed 
from Eq. 2 and uncertainties estimated to be ±2.5 kC. 



Instrument Phase 



(B) (kC) 
Fe II Cr II 



ESPaDOnS 


0.054 


19.4 


17.2 


ESPaDOnS 


0.111 


21.1 


22.9 


ESPaDOnS 


0.640 


10.7 


11.9 


ESPaDOnS 


0.774 


10.9 


12.8 


HARPSpol 


0.787 


11.4 


12.3 


ESPaDOnS 


0.878 


16.9 


14.8 



factor where the L and S subscripts refer to the lines with 
large and small Lande factors respectively. Preston then 
finds the relationship between the surface field and K to 
be 



Bs{kG) = 0.5 + 7.9K. 



(5) 



The challenge with HD 133880 was to find a pair of spectral 
lines that were not strongly blended with other lines from 
which to measure line widths. Table |4] presents two sets of 
lines for Cr ll and Fe ll from which line widths could be 
measured. 



Although good candidates, these spectral lines are still quite 
blended with adjacent lines in the spectrum. To combat this 
problem, we measured line widths at the continuum by re- 
peatedly fitting Gaussians to the spectral line using IRAF's 
splot function. The average width for each feature is taken in 
our calculation. We were also limited by resolving power and 
SNR in our FEROS spectra. Specifically, the combination 
of the lower resolution of the FEROS instrument and lower 
SNRs in the data meant that we were unable to unambigu- 
ously measure the spectral line widths. As a result, only the 
ESPaDOnS and HARPSpol spectra could be utilised. Re- 
sults are summarised in Table [5] and illustrated in the top 
panel of Figure |4l 



It is encouraging that the magnetic field values measured 
from both Fe ll and Cr ll agree well with one another, sug- 
gesting the variations are indeed real. Due to the high blend- 
ing in the spectral lines and the ambiguity in measuring the 
line widths, we estimate the uncertainties to be ± '^ 2500 G, 
though this value may be optimistic. Although this method 
is not as precise as measuring the splitting of spectral lines 
due to the Zeeman effect, it does provide meaningful con- 
straints on the magnetic field. It is clear from these mea- 
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surements that the surface field varies from about 20 kG at 
the negative magnetic pole to 10 kG near the positive mag- 
netic pole (see below). The derived field is consistent with 
a non-sinusoidal curve (though no definitive conclusions can 
be made given the limited number of data points and size 
of the error bars), much in the same manner to the line-of- 
sight magnetic field, but in the opposite sense: the maximum 
of the longitudinal field curve corresponds to the minimum 
of the surface integrated field curve. More measurements of 
(B) are required to constrain the field further which requires 
more high SNR observations of HD 133880. 



6.2 Magnetic field geometry 



The FLDSRCH.F programme was us ed to obtain the mag- 
netic field geometry for HD 133880 (jLandstreet fc Mathvj 
[20001). This programme has longitudinal and surface mag- 
netic field strengths at four different phases as input, which 
are then used to iteratively search for the best fit i, /3, 
Bd, Bq, and Boct- For a more complete descript ion refer to 
iBailey et al.l ()2011f ). Measurements of (Bz) from iLandstreed 
|[1990|), as well as the four new measurements from this work, 
provide good constraints on the observed variations. We also 
have important limits on {B) that appear to range from 
10 kG to 20 kG near the positive and negative magnetic 
poles respectively (see above). This information allowed us 
to derive a new field geometry: i=55° ± 10°, /3=78° ± 10°, 
Ba = -9600 ± 1000 G, B, = -23200 ± 1000 G, and 
Boct = 1900 ± 1000 G. Uncertainties were estimated by ob- 
serving the deviations from the derived field by changing the 
input parameters of (Bz) and (B) within their uncertain- 
ties. This field corresponds well with our calculated value 
of i = 63° ± 18° and has the added advanta ge of agree- 
ing w ith the alternate geometry proposed by [Landstreeg 
()l990l ) (see above). An illustration of our adopted geom- 
etry is presented in Figure [5] The bottom panel of Fig- 
ure [4] shows the mo del (B^) cu rve over-plotted with (Bz) 
measurements from ILandstreed (1990) and this paper. All 
the measurements agree well with the computed model and 
the measurements ma de using the metallic spectrum agree 
remarkably well with ILandstreed (|l990f l who made mea- 
surements from H/3. The field varies non-sinusoidally from 
about -4.5 kG at phase 0.037 to 2 kG at phase 0.537. The 
field curve is fiat and extended around the positive mag- 
netic pole with a rapid decline towards the negative mag- 
netic pole. The quadrupole-to-dipole ratio, q, of the adopted 
model is 2.4, which clearly indicates the presence of a strong 
quadrupol e componen t . We therefore concur with the con- 
clusion of ILandstreed (Il990l ) that the field of HD 133880 
departs in a significant way from a simple centred dipole. 



7 ABUNDANCE MODEL 

We utilise the FORTRAN pro g ramm e ZEEMAN 
(ILandstreed 1 19881 : iLandstreet et al.l Il989l : IWade et al.] 
[2001) to perform abundance analysis of HD 133880. 
ZEEMAN computes emergent spectra of all four Stokes 
parameters for stars with a permeating magnetic field. It 
assumes a low-order multipole geometry for the magnetic 



field based on the above derived values for i, /3, B^, B,, and 
Boct. Up to ten phases of observed Stokes I spectra can be 
fit simultaneously with the abundance of one element varied 
at a time. For each iteration, ZEEMAN automatically 
adjusts V sin i and the radial velocity of the star to optimum 
values. Currently, up to six rings of uniform abundance on 
the surface of the star can be specified that have equal spans 
in magnetic co-latitude. ZEEMAN produces output for 
the abundance distribution of a given element at different 
magnetic latitudes. In effect, this model provides a rough 
ID map of the variation of abundance from one magnetic 
pole to the other. 

Initially, uniform mean abundances for each phase of 
HD 133880 were found, however, we quickly discovered that 
a multi-ring mo d el wou ld provide a better fit. As explained 
by iBailev et all (120111 ). ZEEMAN was designed to quan- 
tify hemispherical abun dance variations f or stars similar to 
53 Cam (=HD 65339) (|Landstreedll988D . HD 133880 ex- 
hibits large-scale abundance variations over the stellar sur- 
face and, as is evident from the longitudinal magnetic field 
curve, both magnetic hemispheres are observed. Our pre- 
liminary abundance analysis suggested that a maximum of 
three rings should be used to quantify the abundance vari- 
ations because we saw abundance anomalies that varied be- 
tween both magnetic poles and the magnetic equator. Nev- 
ertheless, we experimented with fewer and more rings on 
the surface of the star. We found that more than three rings 
did not significantly improve the quality of fits and fewer 
rings did not adequately model the observed variations in 
the Stokes / spectra. 

For our analysis, Tcff was set to 13000 K and log g to 
4.3 which are both consistent with the values determined in 
Sect. 4. The ato mic data were taken f r om Vienna Atomic 
Line Database dPiskunov et al.l 1 19951 : iRvabchikova et al.l 
I1997I : iKupka et al.l Il999l . bOOof T The best-fitting value of 
vsini was found to be 103 ± 10 kms~^. Due to the strong 
blending in the spectral lines, identifying the best fit re- 
quired visual evaluation of the agreement between the com- 
puted and the observed Stokes / spectral lines in several 
windows of the order of 80 A wide. Since the strong mag- 
netic field should suppress any convective motions, which in 
any case are expected to be small at this T^s, the microtur- 
bulence parameter was set to kms^^. 



7.1 Spectra normalisation 

For stars in which vsiui is greater than around 100 kms~^ 
(as is the case for HD 133880) the dominant source of er- 
ror in abundance analysis is from continuum normalisation. 
The reason for this is that in a spectral region that is on the 
order of 100 A wide, there may only exist a few continuum 
points from which to normalise the spectrum. Often it is dif- 
ficult to determine which features should be placed near the 
continuum and which should not. To combat this problem, 
we performed continuum normalisation of spectra of nor- 
mal, sharp-lined, A-type stars from both the ESPaDOnS 
and FEROS instruments to determine the highest degree 
polynomial required to normalise a given spectral window. 
All spectral windows were about 100 A wide to maximize the 
number of continuum points present in the broad-line spec- 
tra. In all cases, no greater than two segments of a cubic 
spline were required to normalise the normal A-star spec- 
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Figure 5. The geometry adopted for HD 133880. The vertical axis is the rotation axis of the star. The angle between the line of sight 
and the rotation axis is j = 55°. The angle between the rotation axis and the magnetic field axis is /3 = 78°. The left panel depicts an 
orientation where the positive magnetic pole is closest to the line of sight (0 = 0.537) and the right panel is when the negative magnetic 
pole is nearer alignment with the line of sight {<j> = 0.037). The bands used for dividing the abundance distribution are shown as solid 
lines and are normal to the magnetic field axis. 



tra correctly. For a given instrument and spectral window, 
we also used two segments of a cubic spline to continuum 
fit the spectra for HD 133880. In this manner, we reduce 
the uncertainty introduced in the placement of the contin- 
uum and ensure that all spectra are normalised in the same 
way. Where necessary, we iteratively performed continuum 
fitting for spectral windows for which synthetic spectra from 
ZEEM AN agreed poorly with our initial normalisation. Hill 
l|l995r) discusses further challenges in performing abundance 
analysis on broad-lined stars. 



7.2 Choice of magnetic field model 

As described in Section 6, two magnetic field models ade- 
quately characterise the observe d variations of the longitu- 
dinal magnetic field, {B^): that of lLandstreetl (|l990l ) and the 
one presented in this paper. To explore the effect of differ- 
ent magnetic field geometries on the derived abundances, we 
performed abundance analysis of HD 133880 for all elements 
using both geometries. The maximum abundance differences 
between the two models were found to be 0.1-0.2 dex. We 
therefore adopt the new magnetic field geometry presented 
in this paper because it agrees well with the value of i derived 
from physical parameters and has the added constraints of 
the observed surface magnetic field modulus variation, (B). 
We report all abundances using the new magnetic field ge- 
ometry. 



8 ABUNDANCE ANALYSIS 



A total of 12 spectra, well spread in phase, were available for 
HD 133880, from the FEROS, ESPaDOnS and HARPS in- 
struments. The HARPS spectrum was not used in the abun- 
dance analysis because its phase corresponded closely to two 
other spectra with comparable SNR and thus added no new 
relevant information. The large wavelength coverage of these 
instruments has allowed us to derive abundance distribu- 
tions for O, Mg, Si, Ti, Cr, Fe, and Pr, as well as upper limits 
for the abundances of He, Ni, and Nd. The large vsini pre- 
sented challenges in deriving abundances due to the presence 
of significant blending in the spectra. Nevertheless, multiple 
relatively clean lines of the majority of elements studied were 
found with varying strengths and Lande factors, providing 
better constraints on the derived abundances. 

As discussed in the previous section, the present dataset 
is sufficient to perform abundance analysis using a three-ring 
model with rings encompassing both magnetic poles and the 
magnetic equator with equal span in co-latitude (60°). Ini- 
tially six spectra were used (2 ESPaDOnS and 4 FEROS) 
that were well spaced in phase to derive the three-ring abun- 
dance model. Once the mean abundance values were found 
for each ring, model fits were produced for all spectra. We 
have no observations for HD 133880 very close to the posi- 
tive magnetic pole (phase 0.537). We do, however, have two 
spectra sufficiently near the positive pole, at phases 0.438 
and 0.640, adequate for providing a first approximation to 
the abundance distributions near the positive pole. 
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Table 6. Abundance distribution of elements studied. Only upper limits of He, Ni and Nd are found. 



He 



Mg 



Si 



Log(nx/nH) 
Ti Cr 



Fe 



Ni 



Pr 



Nd 



Ring — 60° (negative polo) 

Ring 60 — 120° (magnetic equator) 

Ring 120 - 180° (positive pole) 



-1.90 



-2.78 


-3.84 


-2.79 


-5.39 


-4.55 


-3.32 


^ -4.40 


-6.52 


-6.55 


-3.38 


-4.37 


-2.64 


-5.62 


-4.54 


-3.44 


- 


-6.72 


-6.47 


-3.10 


-3.91 


-3.93 


-6.75 


-5.20 


-4.11 


- 


-7.28 


-8.22 


±0.3 


±0.3 


±0.2 


±0.2 


±0.15 


±0.1 


- 


±0.2 


±0.3 



Solar abundance 



-1.07 



-3.31 



-4.40 



-4.49 -7.05 



-6.36 



-4.50 



-5.78 



-11.28 



-10.58 



# of lines modelled 



The mean abundances derived for each element in each 
ring are tab ulated in Table P alo ng with solar abundances 
reported by I Asplund et al.l l|2003 ). The quality of fits were 
tested (where possible) using multiple spectral windows with 
derived uncertainties estimated from the observed window 
to window scatter of the abundance. For lines for which there 
was only one spectral line or region, we estimate the uncer- 
tainties to be about ±0.3 dex, which is the average deviation 
from the derived abundance required to make the quality of 
fit unsatisfactory between the model and observed spectra 
upon visual examination. We also derived a set of abun- 
dance s using Teff = 12000 K as reported by iNetopil et al.l 
([2003) for each of the elements listed in Table [G] We found 
that the abundances in each of the three rings changed by 
±0.10 dex, or less, for the Fe-peak elements (Ti, Cr, Fe) and 
no more than ±0.15 dex for the rare-earths (Pr and Nd), O, 
Mg, and Si which are within our quoted uncertainties. This 
suggests that the choice of T^s does not significantly influ- 
ence the derived abundances for HD 133880. To illustrate 
the quality of fits of the three-ring model, two 80 A win- 
dows for all 11 modelled spectra are shown in Figures [S] and 
[7] The quality and consistency of the model are discussed 
below for each element modelled. 



8.1 Helium 

Many possible lines of He I are considered when determin- 
ing the abundance of helium: 4471, 5015, 5047, 5876, and 
6678 A. However, none of these lines are unambiguously de- 
tected. Nevertheless, fitting these regions provides a useful 
upper limit to the abundance of helium. Our analysis con- 
firms that HD 133880 is a He- weak star with the abundance 
of helium at least a factor of 8 below the solar abundance. 



8.2 Oxygen 

HD 133880 is classified as a Bp star with the A4200-Si 
anomaly which suggests that the oxygen abundance should 
be near or less than the solar abundance (|Sargent fc Searld 
1 19621 : iRobv fc Lambert! llOgd ). Initially, we attempted to 
derive the oxygen abundance using the 7771-75 A triplet. 
These lines suffer from strong non-LTE effects and are highly 
saturated, making them unreliable for abundance determi- 
nation. A more useful set of lines are those at 6155-56-58 A; 
the abundance distribution reported in Table [6] are from 
these lines. These lines are unambiguously detected, but are 



somewhat blended with Pr ll near the redward wing. The 
O abundance varies from the magnetic poles to equator by 
a factor of about 2-3. At all phases, oxygen is found from 
these lines to be overabundant near the negative and posi- 
tive magnetic poles by 0.5 dex and 0.3 dex compared to the 
solar ratio, respectively. Near the magnetic equator, oxygen 
is comparable in abundance to the Sun, within uncertainties. 

8.3 Magnesium 

Only one clean line of sufficient strength to model was found 
for Mg 11 at 4481 A. This line shows some variability, im- 
plying an overabundance of around 0.5 dex near both mag- 
netic poles and a roughly solar abundance near the magnetic 
equator. 



8.4 Silicon 

A total of four lines of Si ii were found that were suitable for 
modelling: 5041 A, the doublet at 5055-56 A, 5955 A and 
5978 A. The first two lines of Si ii were fit simultaneously and 
the result tested by using the same derived abundance model 
in the latter two lines. Satisfactory fits were obtained at each 
phase for all lines of Si, however the doublet at 5055-56 A 
was systematically too strong in the model at all phases (see 
Figure [SJ . Finding the best fit model for the longer wave- 
length lines did not rectify the problem of stronger spectral 
lines in the model and produced a similar abundance distri- 
bution. We find large-scale abundance variations for Si on 
the surface of HD 133880 from pole to pole. At the negative 
magnetic pole and magnetic equator. Si is approximately 25 
times more abundant than in the Sun. The lowest abundance 
is found at the positive magnetic pole where, however. Si is 
still about 3 times more abundant than in the Sun. 

Figure [7] highlights three Si ill lines at 4552, 4567, and 
4574 A that are not modelled well using the abundance de- 
rived from the Si ll lines discussed above (red-dashed line). 
An enhancement of the Si ii abundance of about 1 dex is 
necessary to model well the Si ill lines (blue-dashed line). 
This discordance is very suggestive. The effects of non-LTE 
in the Si ill line s are expected to be sm all at T^s less than 
about 15000 K. iBecker fc ButleJ l|l990l ) demonstrate that 
non-LTE effects in the Si ill lines decrease towards lower 
temperatures. The most probable cause is strong stratifica- 
tion of Si in the atmosphere of HD 133880 with the abun- 
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Figure 6. Spectrum synthesis of tlie region 5000 - 5080 A using a three-ring model. The observed spectra are in black (solid) and the 
model fits are in red (dashed). From top to bottom phases 0.054, 0.111 (both ESPaDOnS), 0.139, 0.325, 0.438 (all FEROS), 0.640, 0.774 
(both ESPaDOnS), 0.809 (FEROS), 0.878 (ESPaDOnS), 0.917, and 0.990 (both FEROS). The line of sight is closest to the negative 
magnetic pole at phase 0.037 and the positive magnetic pole at phase 0.537. Spectra are positioned at arbitrary positions along the 
vertical axis for display purposes. 



dance high in the atmosphere being much lower than near 
an optical depth {tc) of about 1. However, the large usinz 
(103 kms~^) makes normalisation difficult (see Sect. 7.1) 
and studies on more hot Ap stars with low usini's are nec- 
essary to determine how widespread this Si ll/lll discrepancy 
is, and how it may be explained. 



8.5 Titanium 

We have modelled the abundance of Ti using lines of Ti ll 
at 4563 and 4571 A as well as at 4805 A. The final model 
was determined using the former two lines and tested using 
the latter line. Ti exhibits drastic variations between mag- 
netic poles of nearly 1.5 dex which can be seen in Figure [7] 
where lines due to Ti ll are much weaker at phases 0.423 
and 0.625 than near phase 0. Ti is most abundant at the 
negative magnetic pole where it is more than 40 times more 
abundant than in the Sun and transitions to slightly lower 
abundances at the magnetic equator (~ 25 times overabun- 
dant). At the positive magnetic pole Ti is only about 2 times 
more abundant than the solar value. 



8.6 Chromium 

Several lines of Cr were found that were suitable for mod- 
elling, particularly the Cr ll lines at 4558, 4588, and 4592 A. 
The latter two of these three lines were relatively clean while 
A4558 is blended with Fe ii. All three lines were used to com- 
pute the three-ring model presented in Table [B] All Cr lines 
in Figure [7] are fit well at all phases with the model fit at 
4592 A being slightly too strong at later phases. Large vari- 
ations are observed from hemisphere to hemisphere. Cr is 
weakest at the positive pole, but its abundance is still 10 
times greater than the solar abundance. The negative pole 
and magnetic equator have a mean abundance that is ap- 
proximately 100 times greater than in the Sun. 



8.7 Iron 

A plethora of lines are available from which to derive the 
abundance distribution of Fe. The model was found using 
Fe II lines at 4541, 4555, and 4583 A. The lines were of 
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Figure 7. Same as Figure 6 for tlie region 4530 - 4600 A. The green boxes highlight the three Si III lines at 4552, 4567, and 4574 A. The 
red-dashed lines indicate the model fits to abundances derived from the Si II lines that are presented in Table |6] The blue-dashed lines 
are model fits with about a 1 dex enhancement to the Si abundances shown in Table [6] 



sufficient strength for abundance analysis and only the line 
at 4555 A is blended with Cr ll. Similar to Cr, Fe is least 
abundant at the positive magnetic pole, where an abundance 
of about 2.5 times higher than the solar abundance is found, 
and most abundant near the negative magnetic pole and 
magnetic equator where a mean abundance approximately 
10 times larger than in the Sun is observed. Most lines of Fe 
are modelled well, as seen in Figure[7] Some weaker lines are 
systematically too weak in the model (see Figure |6]) which 
is most likely attributable to unrecognized blending. 



8.9 Praseodymium 

Three lines of Pr were found that were useful for modelling. 
The final model was computed using Pr ill and Pr ll at 
6160 and 6161 A respectively. These lines are blended with 
the redward wing of the O I triplet at 6155-56-58 A, but 
are of sufficient strength to obtain a useful abundance. The 
model fits were then tested using Pr ill at 7781 A. A mean 
abundance was found for the two rings covering the negative 
magnetic pole and magnetic equator that is about 10^ times 
overabundant compared to the solar ratio. At the positive 
magnetic pole, Pr is least abundant but still 10* times more 
abundant than in the Sun. 



8.8 Nickel 

There are a few possibly useful lines of Ni ll to model at 5058, 
5059, and 5064 A, however none are unambiguously de- 
tected. Nevertheless, the region (labelled in Figure|6]) proved 
useful for obtaining an upper limit to the abundance of Ni, 
which suggests that it is less than 10 times more abundant 
than in the Sun. 



8.10 Neodymium 

A possible useful line of Nd ill is at 6145 A, however the SNR 
is such that it is not unambiguously detected at all phases 
and rotational broadening causes this line to be strongly 
blended with Fe ll at 6147 and 6149 A. We therefore turn to 
Nd III at 4711-12-14 A. These lines are relatively clean, pro- 
viding useful abundance determinations that reveal a strong 
variation in the mean abundance of Nd between magnetic 
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hemispheres. For both the negative magnetic pole and mag- 
netic equator Nd is greater than 10"* times more abundant 
than in the Sun, whereas at the positive pole this rare-earth 
is about 200 times the solar ratio. 



9 MAGNETOSPHERE 



iLim. Drake fc Linskvl (|l996l ) demonstrated that the 3.5 cm 
and 6 cm radio flux and circular polarisation of HD 133880 
vary significantly and coherently according to the ~ 0.877 d 
period. They reported that the emission shows broad peaks 
near the phases of the longitudinal field extrema (corre- 
sponding, in Landstreet's model, to the phases at which 
the poles of the dipole contribution to the field curve pass 
closest to the line-of-sight), and narrower peaks at the pre- 
dicted phases of quadrupole component pole passages (i.e. 
the quadrupole contribution to the field curve). They also 
reported that the strong extraordinary-mode circular po- 
larisation varied in sign with phase, in agreement with the 
sign of the longitudinal field. As described in Sect. 2, we 
have re-reduced the ATCA radio data, but there are no sig- 
nificant diffe rences in the reduced data compared to that 
described by iLim. Drake fc Linskvl (119961 ). Here we exam- 
ine the radio emission in the context of a magnetospheric 
picture co mbinin g the dynamic a l conc epts summarised by 
lOwocki. Townsend fc Ud-Doulal 11200611 and the radi o mag- 
netosphere model of iLinskv. Drake fc BastianI (|l992r ). 

Using the ste llar wind parameters determine d using the 
CAK formahsm l|Castor. Abbott fc KleinI l|l97d ): values of 
M = lO"^'^ Mq/jt and Voc = 750 kms~^), we compute the 
wind magnetic confinement p arameter ?7« = Bl„ B? / Mvqo — 
10^ (neglecting clumping, lud-Doula fc Owockil l|2002l )'). 
where _Beq is the equatorial surface strength of the magnetic 



dip ole component. The rotation parameter W ■ 



j/Wcrit = 



0.3 l|Ud-Doula. Owocki fc Townseiidl (|2008h . where we have 
included a correction for the oblateness of the star due to 
its rapid rotation). This places the Alfven radius at Raw = 
r}J 7?« ~ 60 -R*. In contrast, the Kepler (or corotation) ra- 
dius is located relatively close to the star, at about -Rkgp = 
VK"^/^ R, = 2.2 R,. This geometry, in which i?Kep << Ram, 
results in a large spatial volume in which magnetospheric 
plasma is predicted to be maintained in rigid rotation with 
the star by the magnetic field, while simultaneously being 
dynamically s upported against gravitational infall d ue to 
rapid rotation (|Ud-Doula. Owocki fc Townsendll2008l ). 

In principle, this situation provides conditions suitable 
for significant accumulation of rigidly-rotating stellar wind 
plasma (diagn osed from optical and UV line emission and 
variability; e.g. lPetit et al.l l|201 j l). as well as a large region 
over which electrons can be accele rated (thus producing the 
non-th ermal radio emission; e.g. iLinskv. Drake fc BastianI 
|l993)). 

The phased radio light curves, illustrated in Fig. [§1 ex- 
hibit a number of interesting characteristics. First, it is clear 
that the radio measurements in Stokes / and V/I vary sig- 
nificantly, and are coherently phased with the rotational pe- 
riod derived from the photometric and magnetic measure- 
ments. Comparison of the 3.5 cm and 6 cm polarisation 
and light curves fully confirm the details of the variations. 
As described by Lim et al., the flux variation is complex, 
characterised by strong, broad maxima at phases 0.0 and 




Figure 8. Lower frame: ATCA Stokes / flux variation at 6 
cm (filled symbols) and 3 cm (open symbols). Upper frame: 
Stokes V/I polarisation at 6 cm (filled symbols) and 3 cm (open 
symbols). 



0.5 (i.e. the extrema of the longitudinal field), and sharper, 
somewhat weaker secondary extrema at quadrature phases 
(i.e. 0.25 and 0.75). The polarisation degree varies approx- 
imately sinusoidally, with extrema of ±16% and a mean of 
zero. The new rotational period derived in Sect. 3 brings 
the phases of the radio fiux and polarisation extrema into 
good agreement with the longitudinal field and photomet- 
ric extrema, resolving the ~ 0.05-cycle offset pointed out by 
iLim. Drake fc Linskvl (|l996l ). The approximate reflectional 
symmetry of the light curve in both I and V/I (horizontally, 
relative to phase 0.0, as well as vertically for Stokes V/I, rel- 
ative to the V/I = 0) is in good qualitative agreement with 
the derived stellar geometry, in particular the large value of 
the magnetic obliquity /3 ~ 90°, implying that we view ge- 
ometrically similar regions of both (magnetic) hemispheres 
during a stellar rotation (see also Fig. 5). 

Comparing the variations at 3.5 and 6 cm, we find that 
the amplitudes are nearly identical in both / and V/I. The 
Stokes I variations may display some small differences in 
shape, in particular with the central (phase 0.5) peak at 3.5 
cm being somewhat shallower and broader than at 6 cm. 

The presence of structures in the radio lightcurve at- 
tributable to the quadrupolar component of the field is very 
interesting. While the Alfven radius for the dipole compo- 
nent fall off as r]J , for the quadrupole the fall-off is as -qj 
(Ud Doula, Owocki and Townsend 2008). The "quadrupo- 
lar Alfven radius" for HD 133880, computed assuming the 



magnetic model derived in Sect. 6.2, is -RAif,q — 21 _R*. This 
implies that the radio emission at both 3.5 and 6 cm forms 
significantly closer than 20 R* from the star's surface, if the 
origin of the radio emission is charged particles in the out- 
flowing stellar wind. 

The uhraviolet C iv AA1548, 1550 and Si iv AA1394, 
1402 UV resonance doublets are both sensitive indicators 
of the stellar wind and its structuring and modulation by 
the magnetic fleld jSchnerr fc Henrichal2008l , e.g.). A search 
of the MAST archive indicates that HD 133880 was unfor- 
tunately not observed by the lUE satellite, nor with any 
other instrument with coverage in this spectral region. On 
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the other hand, our optical spectra contain Ha as well as a 
number of Paschen series H lines. An examination of Ha re- 
veals weak variability in the inner ±500 km/s of the profile. 

In Fig. [5] we illustrate the line profile variations that 
are observed in Ha. Fig. [10] shows the phased equivalent 
width of Ha, along with that measured from the Pa 8590 A 
line. A telluric line removal algorithm was applied to each 
profile before analysis to reduce the variations caused by 
such features. While particular care has been taken to keep 
the normalisation of the ESPaDOnS and FEROS spectra 
as consistent as possible, it is difficult to judge whether the 
small fluctuations that are found in the far wings are real or 
an artefact of the normalisation. At phases 0.0 to 0.35 we 
flnd that the line cores show excess absorption relative to the 
average profiles, reaching maximum absorption at phase 0.0. 
Between phases 0.35 to 0.8, the core region of the Ha profile 
shows emission relative to the average profile. Unfortunately, 
our phase sampling is not sufficient to distinguish if these 
features propagate through the profile, as would be expected 
from rigidly rotating magnetically confined plasma. 

The Ha profile variability of HD 133880 might be a con- 
sequence of modifications of the atmospheric structure due 
to the strong chemical non-uniformities we have determined 
to exist in its atmosphere. However, we know that the disc- 
integrated photospheric He abundance is constrained to be 
below 10% that of the Sun, i.e. less than 1% of hydrogen 
by number, at all phases. So even a significant variation of 
the He abundance could have no important effect on the 
atmospheric structure. 

For the He star a Cen (which has a huge variation of 
He abundance over its surface and s trongly variable metal- 
lic lines: Bohlcndcr. Rice fc Hechleij (|201Q ')) Ha varies only 
in the core of the line (over a region perhaps 2 or 3 times 
the width of a typical metallic photospheric line) and the 
variations quite closely mimic the variations observed in the 
star's iron lines, i.e. the core of Ha gets stronger at the same 
time as Fe ll lines are strong. This may be what we are see- 
ing in HD 133880 as well, since the surface field and metal 
lines are weakest when Ha is also weakest and most of the 
variations occur within ±200 km/s of line centre. Given the 
dominance of the quadrupolar field component, it is possi- 
ble that the contribution of magnetospheric material to the 
overall variability of Ha may be reduced somewhat, since 
material will be trapped in a more complex circumstellar 
geometry: likely four or more "clouds" rather than the two 
clouds seen in stars with dominant dipolar fields. For ex- 
ample, the Ha variability of HD 37776 (Grunhut et al., 
in preparation), a B-type star with a complex quadrupo- 
lar field, is quite modest compared to other stars with more 
dipolar fields ( ((Oksala et al. 2011i , e.g. a Ori E;)) despite the 
fact that it has similar physical properties to the other stars. 
Therefore, at present we simply conclude that the origin of 
the H line variations is unclear. 
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Figure 9. Top: Observed Ha line profiles with the mean profile 
indicated as dashed red. Middle: Shown are the continuum nor- 
malised spectra after subtracting the average of all the Ha line 
profiles, displayed in such a way that the profiles are plotted at 
a height that corresponds to the phase of the observation. Also 
included is a dotted line corresponding to a zero difference to high- 
light the line profile variations. Bottom: Shown is the variance 
of the continuum normalised spectra from the mean profile. 



This paper presents our continuing work in the study of 
chemical abundance evolution with time of magnetic Ap/Bp 
stars in open clusters. The goal of this project is to model 
the magnetic field structure, to derive an estimate of the sur- 
face abundance distribution of a number of elements and to 
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Figure 10. Phased equivalent widths measurements of the core 
region of Ha (top) and the total Paschen line at 8590 A from the 
ESPaDOnS (black squares) and FEROS (red triangles) spectra. 



provide a preliminary analysis to serve as a foundation for 
further more detailed modelling. In particular, more spec- 
tropolarimetric observations will allow for a more sophis- 
ticated interpretation of the surface abundance variations 
and m agnetic field structure vi a Magnetic Doppler Imaging 
fMDI: iHskunov fc Kochukhovl |2002l ). 

HD 133880 is a very young, rapidly rotating Bp star 
{vsvai ~ 103 kms~^) that hosts a strong magnetic field. 
The magnetic field variations were modelled using a co- 
linear axisymmetric multipole expansion with i=55° ± 10°, 
/3=78°±10°, Bd = -9600 ±1000 G, B, = -23200 ±1000 G, 
and Boct ~ 1900 ± 1000 G. We used Hipparcos photometric 
measurements together with our own (Bz) measurements, 
as well as photom etric and radio emission da t a pres ented by 
IWaelkenj j 19851 ) and iLim. Drake fc Linskvl (|l996l ) respec- 
tively, to refine the period to P = 0.877476 ±0.000009 days. 
The abundance distributions of He, O, Mg, Si, Ti, Cr, Fe, 
Ni, Pr, and Nd were modelled using three co-axial rings 
around the two magnetic poles and the magnetic equator. 
The Fe-peak elements are overabundant compared to the so- 
lar ratios with strong variations on a global scale observed 
between hemispheres. 

We have twelve Stokes I spectra of HD 133880 spread 
well in phase that sample both magnetic poles and the 
magnetic equator. Our adopted magnetic field geometry, al- 
though an approximation to the actual field structure of 
HD 133880, and our three co-axial ringed abundance model 
provides the most detailed abundance analysis possible given 
the limitations of the model. Most elements studied appear 
to have a distinctly non-solar abundance in most rings. The 
abundance of Mg is slightly above solar at both magnetic 
poles but near solar abundance at the magnetic equator. 
Upper limits to the abundances of He and Ni were found. 
HD 133880 is clearly a He-weak star with an abundance 
at least a factor of 8 lower than the solar ratio. Ni is no 
more than 10 times overabundant compared to the Sun. All 
other elements studied (O, Si, Ti, Cr, Fe, Pr, Nd) are more 
abundant at the negative magnetic pole than the positive 
magnetic pole. Most interesting are the variations observed 
in O abundance. Most Ap/Bp stars s how O abundances near 
or below the solar abundance ratios (jSargent fc Searlell 19621 : 



iRobv fc Lamberjl 199(3 ). However, a careful study of several 
O features reveal that oxygen is more abundant than so- 
lar at both magnetic poles (most notably at the negative 
magnetic pole where O is 0.7 dex overabundant) but a so- 
lar abundance near the magnetic equator. This behaviour 
is certainly uncommon in Bp stars, but is unambiguously 
determined in this analysis. 

In the context of the larger study of the time evolution 
of atmospheric abundances of Ap/Bp stars, it is interesting 
to compare to work done bv lBailev et al.l (|201ll ) on the Ap 
star HD 318107 (=NGC 6405 77). This star is also hot {T^g 
~ 11800 K), young (logi ~ 7.8 (yrs)), of comparable mass 
(2.95 ± 0.15 Mq), and host to a large complex magnetic 
field {{Bz) varies from about 1 to 5 kG). In contrast to 
HD 133880, it is a slow rotator {vsini ~ 7 kms~^) and only 
one magnetic hemisphere is observed. The global properties 
of atmospheric abundances is similar between both stars: the 
mean abundances for Fe-peak (Ti, Cr, Fe) and rare-earth 
(Nd, Pr) elements are comparable within uncertainties and 
clearly overabundant compared to the solar ratios. There 
are, however, notable differences in the abundances of O and 
Si between both stars. For HD 318107, the O abundance is 
more typical of Ap stars, being about 0.3 dex lower than 
the solar value as compared to HD 133880 where O appears 
to be overabundant. Si is about 1 dex more abundant in 
HD 133880 compared to HD 318107 which is understandable 
in that the former has the Si A4200 peculiarity whereas the 
latter does not. 

iFolsom et al.l (l2007|) studied the chemical abundances 
of stars in an older cluster NGC 6475 (logt ~ 8.47 (yrs)) 
that have similar masses to HD 133880. Mean abun- 
dances for multiple elements in the Ap stars HD 162576 
(=NGC 6475 55; 3.0 ± 0.4 Mq ), HD 162725 (=HR 6663; 
3.3±0.5 Mq) and HD 162305 (=NGC 6475 14; 2.7±0.3 Mq) 
were derived. As compared to HD 133880 and HD 318107, 
the abundances of Cr for all three stars are similar, being 
about 2 dex larger than in the Sun. However, the abundances 
of Ti and Fe appear to be slightly lower and only 0.5 dex 
larger than the solar ratios, respectively. This is in contrast 
to HD 133880 and HD 318107 where Ti and Fe are about 1.5 
and 1 dex overabundant compared to the Sun, respectively. 
A larger dataset is required before any definitive conclusions 
can be made regarding the time evolution of abundances in 
Ap stars. 

HD 133880 was previously modelled by [Landstreeg 

(|l990l) who characterised the magnetic field variations of this 
star with measurements of the longitudinal magnetic field 
from H/3. Unlike most magnetic Ap/Bp stars whose mag- 
netic field structure is predominantly dipolar, iLandstreeti 
(jl990l ) found a magnetic field dominated by the quadrupole 
component. Despite the large v sin i of the star, we were able 
to further characterise the magnetic field of HD 133880. The 
higher SNR and spectral resolution of the ESPaDOnS spec- 
tra made it possible to compare spectral lines of the same 
element that have a large and small Lande z factor to char- 
acterise the surface magnetic field modulus variations, (B). 
By finding a model that characterises the observed (Bz) and 
(B) variations, we are able to find the inclination, i, of the 
line of sight to the rotation axis as well as the tilt angle of 
the magnetic field axis, /?, to the rotation axis. This model 
is merely an approximation to the very complex field struc- 
ture of HD 133880, but is sufficient to provide a meaningful 
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Table 7. Measurements of (Bz) using LSD masks composed of 
lines of single elements (Si, Ti, Cr, and Fe) from the polarised 
ESPaDOnS and HARPS spectra (cf Table [g. 



Phase 


(B,)(Si) 


(-B.)(Ti) 


(S.XCr) 


(S,){Fc) 




(G) 


(G) 


(G) 


(G) 


0.111 


-3419 ± 76 


-5724 ± 535 


-4773 ± 229 


-3498 ± 50 


0.640 


1790 ± 66 


2255 ± 304 


2653 ± 196 


1993 ± 44 


0.787 


1020 ± 93 


-506 ± 308 


1026 ± 208 


975 ± 78 


0.878 


-1776 ± 63 


-2883 ± 319 


-2409 ±175 


-2010 ± 45 



first analysis to tlie atmosplieric abundance distribution of 
elements. 

Figure Hlsiiows clearly the qualitative nature of botii tiie 
(Bz) and {B) variations. There is a rapid change in polarity 
in the {B^) curve with a steep transition from the positive 
to negative hemisphere. The strong quadrupolar component 
of the magnetic field causes a distinct non-sinusoidal varia- 
tion with the field structure flattening out near the positive 
magnetic pole. The fact that there is a clear pole rever- 
sal indicates that the sum of the inclination and magnetic 
field axis must be greater than 90° : i + 13 > 90° , and that 
the line-of-sight goes well into both magnetic hemispheres. 
The parameters of the adopted magnetic field geometry are 
presented in Sect. 6.2. (Bz) measurements using all metal- 
lic lines (see Tab led agree well with the measurements by 
iLandstreed ([19901) using H/3, suggesting that they sample the 
magnetic field in similar ways (see Fig.|4]). Longitudinal field 
(Bz) measurements of spectral lines of only a single element 
provide clues to the inhomogeneous structure of the mag- 
netic field and the surface distribution of elements. Using 
all four of our polarimetric observations, we have remea- 
sured the longitudinal field {Bz) using spectral lines of Si, 
Ti, Cr, and Fe. The results are listed in Table [T] 

It is clear from Table [7] that all four elements sample 
the field in different ways. The differences in measured field 
strengths can b e attributed to th e complex magnetic field 
(as reported by Landstrect (1990) and discussed in this pa- 
per) and the inhomogeneous surface distribution of elements 
which is a common feature of Bp stars. Near the positive 
magnetic pole (phase 0.537) at phase 0.640, as well as phase 
0.878, the measured (Bz) values for all four elements dif- 
fer significantly. In the negative magnetic hemisphere (near 
phase 0.037) Si and Fe sample the field in similar ways, 
whereas Cr and Ti sample the magnetic field differently with 
{Bz) measurements approximately 1 and 2 kC stronger than 
Si and Fe. Along the sharp decline from the positive to neg- 
ative magnetic pole (near phase 0.787; see Figure [4|, it ap- 
pears that only Ti differs greatly showing a polarity reversal 
as compared to the other elements at the same phase. 

The measurements of the longitudinal field using dif- 
ferent elements reveal that particular elements at different 
phases sample the field in very distinctive ways which sug- 
gests a truly complex magnetic field that cannot be accu- 
rately modelled by a simple axisymmetric ma gnetic field 
(such as the one adopted in this paper and by iLandstreeJ 
|1990|)). It also strongly suggests an inhomogeneous distri- 
bution of various elements on the stellar surface that sub- 
stantially affect the magnetic field measurements. We also 
compared the computed to the observed Stokes V profiles 
for the three polarimetric ESPaDOnS spectra for multiple 



lines (Mg, Si, Ti, Cr, and Fe). The agreement is poor, with 
the computed Stokes V profiles being much stronger than 
the observed, further suggesting that our adopted magnetic 
field model is only a rough a pproximation. Th i s resul t is 
not unique. A recent paper bv lKochukhov et al.l l|2011al ) on 
HD 37776 (=HIP 26742) reveal that the complex multipole 
model derived from the longitudinal field predicts Stokes V 
signatures that are much too strong compared to the ob- 
served Stokes V profiles. In this case, the multipole mag- 
netic field model (similar in structure to the one used for 
HD 133880) predicts a large surface field, whereas the MDI 
results produce a much more modest surface field. This is 
also seen, although less severely, for a Ori E (jOksala et al.l 
I2OIH ). For this star, detailed analysis using MDI produced 
a magnetic field model that closely resembled the observed 
Stokes V profiles. The acquisition of new, high SNR polari- 
metric observations of HD 133880 are necessary to charac- 
terise unambiguously the surface abundance variations and 
magnetic field structure using MDI. 

The very strong field and rapid rotation of HD 133880 
produce interesting magnetospheric characteristics. In our 
preliminary analysis, we were not able to conclude that the 
optical spectra contain detectable information about the 
magnetosphere. However, HD 133880 is one of only a few 
Bp stars for which radio flux variations have been noted. 
There are small but distinct differences between the 3 cm 
and 6 cm radio emission, such as the systematically lower 3 
cm than 6 cm fluxes at phases 0.0 and 0.5. The 3 cm emis- 
sion suggests a "bump" at about phase 0.6 and the peak at 
about phase 0.75 is narrower at 3 cm than at 6 cm. The fact 
that the radio spectrum is flat between 6 and 3 cm suggests 
that the emission is optically thick, and therefore we are 
observing the magnetosphere at two different radio surfaces 
rather than an integration throughout the entire volume. 
These subtle differences in the radio light curves at 6 and 3 
cm then imply that these two surfaces cannot have exactly 
the same shapes and that, presumably, the magnetic fleld 
is not perfectly axisymmetric. HD 133880 is a very exciting 
target for more detailed magnetospheric modelling. 
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